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Рассмотрены обоснование и синтез математической модели линий электропередачи, позволяющей воспроизводить весь спектр 
нормальных и аварийных процессов в трехфазных линиях электропередачи различной протяженности с учетом электромагнит- 
ного взаимовлияния параллельных цепей и возможного коронирования. 


Исходной основой при любом способе матема- 
тического моделирования электроэнергетических 
систем (ЭЭС) служат математические модели эл- 
ементов ЭЭС, от точности которых зависит досто- 
верность результатов моделирования. Нередко дру- 
гих критериев ее оценки, особенно при воспроиз- 
ведении сложных аварийных процессов, нет. 

В сложившейся практике математического мо- 
делирования реальных ЭЭС линии электропереда- 
чи (ЛЭП) моделируются алгебраическими уравне- 
ниями П-образной схемы замещения. 

Значительное усложнение современных ЭЭС и 
насыщение их быстродействующими средствами 
релейной защиты, технологической и противоава- 
рийной автоматики, а также внедрение более бы- 
стродействующего коммутационного оборудова- 
ния, ставит задачу повышения точности моделиро- 
вания аварийных процессов в ЭЭС. Причиной это- 
му служит также постоянно значительное число тя- 
желых аварий, из-за недостаточно полной и досто- 
верной информации моделирования, используе- 
мой при проектировании и эксплуатации. 

Поэтому разработка математической модели 
трехфазных ЛЭП, обеспечивающей адекватное 
воспроизведение непрерывного спектра нормаль- 
ных и аварийных процессов в ней является весьма 
актуальной. 

В отличие от прочих элементов ЭЭС, физиче- 
ская сущность которых позволяет обоснованно 
считать их элементами с сосредоточенными пара- 
метрами, для ЛЭП высокого и сверхвысокого на- 
пряжения длиной /> 50 км целесообразен учет ра- 


спределённости параметров, а также нередко воз- 
никающего коронирования, которое влияет на по- 
тери в ЛЭП и на электромагнитные процессы в це- 
лом. Еще одним фактором, который следует учиты- 
вать в математической модели ЛЭП, является элек- 
тромагнитное взаимовлияние параллельных цепей. 

Система телеграфных уравнений и электрофи- 
зическая схема замещения малого участка или /-ой 
ЛЭП в целом определяют исходной математиче- 
ской основой любой модели ЛЭП, общеизвестную 
систему уравнений: 
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где: Іф, - индуктивность контура: фаза линии - об- 
ратный эквивалентный провод (земля), обусловлен- 
ная только магнитным потоком в воздухе, при пре- 
небрежении, ввиду малости, составляющей индук- 
тивности, определяемой магнитным потоком в са- 
мом проводе фазы, и в предположении бесконечно 
большой проводимости земли; Ь ЗІ - дополнительная 
индуктивность контура: фаза-земля, обусловленная 
только магнитным потоком в земле; М Ші - взаимо- 
индуктивность между любыми двумя из трёх конту- 
ров: фаза-земля; г фі - активное сопротивление про- 
вода фазы линии; г ЗІ - активное сопротивление об- 
ратного эквивалентного провода (земли); С ф , - ём- 
кость провода фазы относительно земли; С МФ; - ём- 
кость между любыми двумя проводами фаз линии; 
і м , і к , і а - мгновенные значения фазных токов; 
и м , и ВІ , и а - мгновенные значения фазных напряже- 
ний; 4, = (/ф+/в,+/'а)Л _ мгновенное значение тока ну- 
левой последовательности; и ОІ =(и АІ +и ВІ +и СІ )/3 - 
мгновенное значение напряжения нулевой последо- 
вательности; /, - длина /-ой ЛЭП или её участка. 

Замена правых частей уравнений системы ко- 
нечными приращениями напряжений и токов по 
длине ЛЭП приводит к её модели в виде конечного 
числа трёхфазных П-образных схем замещения, ко- 
личество которых определяется длиной ЛЭП и тре- 
буемым уровнем адекватности воспроизводимых 
процессов. Так как при таком методе решения каж- 
дый из участков линии отображается аналогичной 
системой уравнений с модифицированными правы- 
ми частями в форме разностей соответствующих на- 
пряжений и токов, то в целом математическая мо- 
дель длинной линии оказывается весьма громозд- 
кой. Поэтому подобный способ моделирования 
оправдан для сравнительно коротких ЛЭП. Прие- 
млемая для практики точность воспроизведения пе- 
реходных процессов указанным способом обеспе- 
чивается при длине участка или линии в целом 
/,<250 км. Более строгое моделирование процессов в 
ЭЭС, в частности связанных с коммутационными 
перенапряжениями, предполагает дальнейшее уме- 
ньшение этого параметра до значений /,<50 км. 

В общем случае, в связи с необходимостью учё- 
та магнитного и электростатического взаимовлия- 
ния параллельных линий или соответствующих 
участков линий, а также грозозащитных тросов, 
которое по данным многочисленных исследований 
проявляется преимущественно в токах и напряже- 
ниях нулевой последовательности, оптимальной 
является запись уравнений в системе несимме- 
тричных составляющих а, Д 0. Достоинством та- 
кой математической модели, наряду с её компакт- 
ностью, является естественное, а не дополнитель- 
ное, наличие уравнений, связывающих параметры, 
напряжения и токи нулевой последовательности. 
Последнее позволяет наиболее просто и удобно 
осуществлять учёт электромагнитного взаимовлия- 
ния параллельных цепей. 

Для некоторой /-ой ЛЭП или её участка исход- 
ная система уравнений в системе а, Д 0 приобрета- 
ет вид [1]: 
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- представляют собой мгновенные значения токов 
и напряжений в системе а, Д 0; Д и С, - индуктив- 
ность и ёмкость для прямой и обратной последова- 
тельностей, именуемые также рабочей индуктив- 
ностью и ёмкостью трёхфазной линии; г, - актив- 
ное сопротивление трёхфазной линии для прямой 
и обратной последовательностей, представляющее 
собой активное сопротивление провода фазы; Ь ш , 
С 0І и г ы - собственные индуктивность, ёмкость и 
активное сопротивление /-ой трёхфазной линии 
для нулевой последовательности. 

При этом мгновенные значения фазных вели- 
чин связанны с составляющими а, Д 0 формулами: 
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Используемые в уравнениях реактивные и ак- 
тивные параметры ЛЭП определяются многократ- 
но апробированными выражениями, которые де- 
тально и с учётом возможных исполнений трёхфаз- 
ных линий рассмотрены в [2]. Особое внимание в 
этих и других аналогичных работах уделяется пара- 
метрам нулевой последовательности в связи со 
сложной спецификой индуктивного и электроста- 
тического взаимодействия линии и земли, а также 
сравнительно близких параллельных цепей: линий, 
участков линий, грозозащитных тросов и т.д. В 
частности, результаты исследований и опыт дли- 
тельной эксплуатации ЛЭП высокого и сверхвысо- 
кого напряжения указывают на целесообразность 
учёта для нулевой последовательности явления по- 
верхностного эффекта в земле, тогда как в прово- 
дах линии им можно пренебрегать. 

Активное сопротивление земли, определяемое 
формулой: г,=10~Ѵ/ Ом/км, практически не зависит 
от удельной проводимости грунта у= КГ 5 . . . 1 (Г 2 Ом/схг 1 , 
т.к. изменение последней компенсируется соответ- 
ствующим изменением сечения, по которому рас- 
пространяется ток в земле. Глубина расположения 
эквивалентного обратного провода, имитирующе - 


207 




Известия Томского политехнического университета. 2005. Т. 308. № 7 


го этот процесс, определяется формулой [2]: 
п 2,085 

О, = =, м, полученной для токов возвра- 

Ѵ/.уіО- 3 

та фиксированных частот, а индуктивность такого 
провода, обусловленная магнитным потоком в зе- 
мле, выражается формулой [2]: 


4, -21п— — -10 4 , Гн/км, 

2 К 


где/ Гц, - частота тока возврата, в выражениях для 
г 3 и Д ; Лер, м, - средняя геометрическая высота про- 
водов фаз над землёй. 

Из формул для г, и Д очевидна непосредственная 
зависимость процесса распространения тока в земле от 
его частоты, а также возможность отображения этого 
процесса с помощью эквивалентного обратного прово- 
да, состоящего из параллельных г 3 Д-ветвей с распреде- 
лёнными или сосредоточенными параметрами. 

Согласно данным [1] моделирование эквива- 
лентного обратного провода тремя параллельными 
/■Д- ветвям и обеспечивает достаточную точность в 
диапазоне частот до 1000 Гц, что охватывает весь 
реальный спектр переходных процессов в ЭЭС, 
включая коммутационные перенапряжения. Рас- 
чёт параметров указанных г, Д-ветвей осуществля- 
ется исходя из условия: 


*= 1 г зк + І со ^тк к л (б) к ) + ісоЬ 3 ((О к ) 


где г л , Ь л - искомые параметры г, Д-ветвей, которые 
определяются в результате решения системы из ше- 
сти алгебраических уравнений (три для г зк и три для 
Ь зк ) для трёх значений со ь обычно: 70, 200, 800 Гц. 

Принимая во внимание уточняющий характер 
учёта поверхностного эффекта в земле и связанное 
с этим усложнение математической модели линии, 
следует считать целесообразным его применение 
только для собственных параметров нулевой по- 
следовательности линии, включая грозозащитные 
тросы, учёт которых не изменяет структуры мате- 
матической модели, а лишь влияет на значения 
указанных параметров. Что касается электромаг- 
нитного взаимовлияния параллельных линий, 
участков, то этот фактор, согласно многочислен- 
ным данным, считается достаточным учитывать 
традиционным способом, предполагающим для 
каждого индуктируемого тока нулевой последова- 
тельности одну г^Ц- ветвь. 

Рассмотренные аспекты математического моде- 
лирования ЛЭП определяют и характеризуют два 
возможных варианта решения этой задачи. В соот- 
ветствии с одним из них основу математической 
модели составляет какая-либо из приведенных вы- 
ше систем уравнений, записанная для конечных 
приращений напряжений и токов в правой части 
уравнений, которая в этом случае отображает про- 
цессы в линии или участке с сосредоточенными па- 
раметрами в виде П-образной схемы замещения. 
При этом, качество моделирования процессов, 


определяемое спектром воспроизводимых гармо- 
ник, зависит, как уже отмечалось, от длины линии 
и связанного с ней количества участков. Посколь- 
ку такая информация априори доступна лишь при- 
менительно к конкретной ЭЭС, то очевидно, что и 
данный вариант, без всегда нежелательной избы- 
точности, приемлем именно для таких условий. 

В общем случае, предполагающем универсаль- 
ность использования математической модели трёх- 
фазной линии безотносительно к её параметрам и 
спектру воспроизводимых процессов, особенно их 
высокочастотной части, наиболее эффективен и 
удобен другой способ применения этих систем ура- 
внений, в частности для а, Р, 0, учитывающий ра- 
спределённость параметров линии электропереда- 
чи и поверхностный эффект в земле. 

Теоретическую основу обозначенного способа 
составляет волновая интерпретация процессов в 
этих цепях, которая позволяет представить исход- 
ную модель системой уравнений: 

(Р) ~ (Р) = [ Ц «О») - 4 Ч к № “ 

Ѵ<ѵ(р) + =[ ЧшЫ + 


где Н и К - обозначение условного начала и конца 
линии или её участка; Е, - индекс составляющих 


системы а, (5, 0; 2,. 



- волновое сопротив- 


ление идеализированной линии для соответствую- 


щих составляющих системы а, /3, 0; Д = (^/44; 

- постоянная изменения фазы волны соответ- 
ствующих составляющих системы а, /?, 0 в идеали- 
зированной линии длиной /,; 


<т с , = 


22, 


- + ^2 а . 
2 * 


- параметр затухания соответствующих соста- 
вляющих системы а, р, 0, обусловленного потеря- 
ми в реальной линии длиной /, из-за гдО и #д0. 

Модификация системы уравнений, необходи- 
мая для учёта поверхностного эффекта в земле, вы- 
текает из содержания выражений, определяющих 
2у, Ту, а у для нулевой последовательности: 2 0/ , т 0/ , 
сг 0 , и рассмотренного ранее способа отображения 
распространения токов нулевой последовательно- 
сти фиксированных частот посредством трёх па- 
раллельных /•„/.„-ветвей: г' 3 Ц, г"Ц, г"1"'. 

Замена в выражениях 2„ г,,, для 2 0 „ г„„ <т 0; , г зі 
и /„ на /•',., Ь приводит к образованию формул для 
параметров нулевой последовательности: 
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и аналогичных для г", Ц, г"\ Ц' и соответственно 
трёх волновых каналов нулевой последовательно- 
сти с этими параметрами. 

Для учета электромагнитного взаимодействия, 
проявляющегося в появлении в напряжениях нуле- 
вой последовательности /-ой параллельной цепи: 

Ѵш (0 = ^,ш (0 + Ѵвш (0 + Цгн .(0]; 

ѵ™ (0 = (0+' (0 + (01, 

кроме составляющих этого напряжения и 

обусловленных распространением собствен- 
ного тока нулевой последовательности /-ой цепи 
Тіш(0 и / ок ,(0, а с учётом поверхностного эффекта в 
земле: / 0 ' ш (/), /он, -(0, /он' (О и 4,(0, %,(/), 4(0, соответ- 
ственно, дополнительных составляющих: 
ад, ^кя(0;-; 0^(0, Ч/К-'МнМ ѴкМ, Образован- 
ных распространением токов нулевой последователь- 
ности: /оні(0,/окі(0;-;гон;(0,гок,(0;-Дош(0,4ю,(0, влия- 
ющих цепей: 1,...,/,...,я по соответствующим взаим- 
ным волновым каналам, математическая модель /-ой 
многоцепной ЛЭП должна содержать, помимо рас- 
смотренных ранее основных уравнений, дополни- 
тельные уравнения, обеспечивающие формирование 
напряжений: Ц т ,({) и Ѵ т {і) указанным выше образом. 

С учетом вышеизложенного, ранее записанная 
система уравнений приобретает вид: 

4н,(0 - ш(0 = [4 Ю 0 -Е)- ^ ю (/ - т,.)] • ; 

4кД0 + г Аш(0 = - 4 + 2,4 н,( ? - т,-)] • 

^ш(0 — зщ(0 = [^/ш0 — Т і) — -^4 К/0 — т ,)] ' е ’ > 

^ю(0 + ^ю(0=[^щ(^-т,-) + ^щ0-т,-)] •е“ СТ; ; 

(0“ 4Л>ш (0 = [^ок, 0 — ^о; ) — ^оАю 0 — 7 о/ )] ' е °' ' > 

4к,- (0 — ^оіЧкі (0 = 14ш0 — т о;) + — т о,)1 ' е 

и далее по аналогии для ^о,, т , , ; , гт м ^, т о,, сх^, в ко- 

торых: 

*оні (0 = Ош (0 + Ош (О + Он 1 (О, 

О к; (0= ОкДО + 0ю (О+Іокі (О, 

4нД0 = ^ш(0 ■ -Е/ И і(0 ■ -••• (0 -■■■ -Ущ, (О; 

4кД0 = і/ ю ,(0-і/ ш (0-... -с/ Ю у (0-- (О; 

4нД0 = I' (0 - ^н, (0 + ДО]; 

Ѵ / т,(» = ^[Ѵш(0-ѴсшЛ0]; 

V шД г ) = -[^ш ( 1 ) + ( 1 ) + 4ш О )]; 

4кД0 = |^к, (0 - з Ѵви (0 + ^ск, (01; 

^кД0=^[^кД0-і/ С к,Д0]; 

^ ю Д0 = ^ ю (0+^ ю (0+і/ С к ,(01; 


1 -\/3 

Ош (0 = 0н> + , аНІ’ ! Ш, (0 — ОнДО “ ^ОнДО + ДО, 

і Дз . 

г снД0 = Он ДО — ~ОнДО ДО; 

1 УІЪ 

ОкДО = *ою +0к,; ОкДО = ОкДО — ^ОкДО + ~^~0к ДО; 
і Дз . 

*скД0 = ОкДО — -ОкДО ДО- 


где 2. = 2 а/ 



- волновое сопротивле- 


ние для составляющих а, Д такое же, как и для со- 
ставляющих прямой и обратной последовательно- 
стей, а/(/-т ; ,) - функция запаздывания. 

Полученная система уравнений учитывает все 
ранее названные факторы, за исключения корони- 
рования ЛЭП, которое связанно с возникновением 
при определенных условиях объемного заряда на 
проводах ЛЭП, описываемого соответствующими 
вольт-кулоновыми характеристиками. Однако, 
если иметь в виду, что коронирование приводит к 
динамическому изменению активной проводимо- 
сти между проводом и землей, и емкости, то влия- 
ние короны можно учесть путем динамической 
коррекции параметров математической модели 
и “С”, в выражениях 


1 


4 

V 

.+ Ёіі-2 

4 г =/ , 

І1 іі с іі и Ч = 


г ’ Т и ч 

22, 

\ ь 

2 


которые вычисляются по формулам [1,3]: 

- динамическая активная проводимость корони- 
рующего провода: 


8 = \ 


Г/Г' 62 

г 

1-ехр 

у50 ) 

V 


-3,05 


V -V 

т н; 

С/ 


10 3 , МОм-м; 


где \] т - мгновенное значение фазного напряже- 
ния; / - частота; 1/ т , - начальное напряжение ко- 
роны, определяемое выражением [3]: 


V =Е Діп— , кВ, 

нач нач п ’ 5 

К 


в котором Н - высота провода над землей, а К- ра- 
диус провода; Е шч — начальная напряженность 
электрического поля при коронном разряде, вычи- 
сляемая по формуле [3]: 


Е иш =24,5тд 


1 + 


0,65 
(<5 К) 0 ’ 3 


кВ/см, 


где м=0,82...0,94 - коэффициент гладкости прово- 
да; 5 - относительная плотность воздуха; 


- динамическая емкость коронирующего провода 

[ 3 ]: 


С = С Г + 2,4 


50 

ѵ/ 


0,42 /■ 


V, 

. V 

\ нач 


— 1 


пФ/м; 


где С г - обычная геометрическая емкость провода. 
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Рисунок. Пример осциллограммы однофазного короткого замыкания в НЭП 500 кВ (ЕЕРЗ_2), где: І-ЕРЗ_2.Іа, І-ЕРЗ_2.ІЬ, 
1ЕРЗ_2.!с - фазные токи; ЕЕРЗ_2.ІЮЗ ~ ток нулевой последовательности 


Заключение 

В результате рассмотренного синтеза получена 
универсальная всережимная математическая мо- 
дель ЛЭП, позволяющая воспроизводить непре- 
рывный спектр нормальных и аварийных процес- 
сов в линиях любой протяженности и различных 
классов напряжения с учетом уточненного элек- 
тромагнитного взаимовлияния параллельных це- 
пей и возможного коронирования ЛЭП. Примене- 
ние трех параллельных волновых каналов для нуле- 
вой последовательности существенно повышает 
точность моделирования наиболее вероятных в эк- 
сплуатации ЭЭС несимметричных аварийных про- 
цессов, включая коммутационные перенапряже- 
ния. 


Разработанная модель может использоваться 
при проведении расчетов как численными метода- 
ми и средствами, так и гибридными. Последнее реа- 
лизовано в гибридных процессорах ЛЭП всережим- 
ного моделирующего комплекса реального времени 
ЭЭС гибридного типа, разработанного в НИЛ «Мо- 
делирование ЭЭС» Электротехнического института 
Томского политехнического университета. 

Экспериментальная проверка гибридных про- 
цессоров ЛЭП в составе моделирующего комплек- 
са Тюменской энергосистемы подтвердила аде- 
кватность синтезированной модели ЛЭП для всего 
спектра воспроизводимых процессов, фрагмент 
которых для наиболее распространенного вида ава- 
рий приведен на рисунке. 
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